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its properties which may be deduced from observations of pulsar glitches and postglitch 
behavior. Current observational evidence concerning the softness or stiffness of the high 
density neutron matter equation of state is reviewed briefly, and the (revolutionary) 
implications of a confirmation of the existence of a 0.5ms pulsar at the core of SN1987A are 
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Introduction, 

Neutron stars are an appropriate topic for a symposium which celebrates the 
remarkable oeuvre and persona of Lev Davidovich Landau, in part, because, according 
to Leon Rosenfeld, it was Landau who first suggested their existence. At a gathering at 
the home of Niels Bohr, on the evening of the day that the news of Chadwick's discovery 
of the neutron had reached Copenhagen, Landau, who had been interested in 
understanding how stars work, proposed that one might have stars made up entirely of 
neutrons. Presumably Landau, who had, independently of Chandrasekhar, recognized 
that the degeneracy pressure of electrons makes it possible for a white dwarf of mass 
less than ~ 1 ,4M© to remain stable against gravitational collapse, had in mind the 
possibility that a neutron star could accomplish the same feat through the degeneracy 
pressure of the neutrons. Some seven years were to pass before Robert Oppenheimer 
and George Volkoff carried out the initial analytic investigation of the stability of neutron 
stars, and, quite independently of Landau, showed that this did indeed occur. 

Additional reasons to discuss neutron stars in a meeting dedicated to Landau 
come from the opportunities which they offer for the study of matter under extreme 
conditions. With average densities in excess of that of nuclear matter (p 0 = 2.8 x 10 14 
gm cm- 3 ), neutron stars are the highest density observable matter. Relatively speaking, 
they contain the lowest temperature superfluid. For example, the neutron liquid core of 
the pulsar, PSR0525-21 , which is some 10 6 y old, is ~ 10 5 K, so that its temperature 
relative to its degeneracy temperature (100 MeV ~ 10 12 K), is lower than that of 
superfluid helium at a few tenths of a microdegree. On the other hand, the neutron 
superfluid which coexists with a periodic lattice of neutron-rich nuclei in the inner crust of 
a neutron star (at densities p such that p o /500 5 p Sp 0 /2), is the ultimate high 
temperature superfluid, since its transition temperature may reach 8 x 10 9 °K. Moreover, 
the neutron superfluids in the crust and core of neutron stars are also by far the most 
abundant superfluids we encounter in the universe. 


Like the work of Landau, the study of neutron stars embraces almost every sub- 
field of physics: low temperature physics, condensed matter physics, nuclear and 
particle physics, gravitational physics, plasma physics, hydrodynamics. The physicist or 
astrophysicist working on neutron stars thus encounters Landau at essentially every 
step in his research. 

On the other hand, his friends and colleagues tell us that Landau disliked 
astrophysics, or at least those branches of astrophysics in which investigators explore 
multiple scenarios, conclude that a large sub-set of these are consistent with 
observation, and present progress reports on these not very satisfying conclusions at 
regular intervals. Indeed, astrophysics differs from physics in that one cannot perform 
experiments to verify theory. Rather the theoretical astrophysicist must attempt to 
develop scenarios for a given body of observational data which are sufficiently detailed 
that future observations can verify or disprove a particular scenario, and provide thereby 
a direct link between theory and observation. In some areas (and these might have 
been of interest of Landau) astrophysics has matured to the point that there exists 
sufficient observational data to rule out large classes of theories and, conversely, to 
persuade us that a given physical picture is correct. In my talk today I shall consider the 
extent to which the study of the interiors of neutron stars has reached this stage of 
maturity, and thus give you a progress report on neutron stars as a cosmic hadron 
physics laboratory. 

A CL Qygr v ifiW 

In Figure 1 I reproduce a cross-section of a 1 .33Me neutron star calculated using 
a Bethe-Johnson equation of state for neutron matter at densities p £ p 0 . The softness 
or stiffness of the equation of state in this region determines the mass-radius relation, 
the crustal extent and the stellar moment of inertia for a given mass, and the central 
density of the star. For this "moderately stiff" equation of state, one has a solid crust 



which extends some 1 .6 km below the surface of the star and a quantum liquid interior 
which contains electrons, superfluid neutrons in a 3 P 2 pairing state, superconducting 
protons in the 1 S 0 pairing state, and possibly, a pion condensate at densities greater 
than 2p 0 . Of particular interest to us will be that region of the inner crust lying between 
densities p o /10 5 p 5 p 0 /2, which current calculations show contains a neutron 
superfluid in the pairing state which is pinned to the lattice of neutron-rich nuclei. 

I summarize some current efforts to use neutron stars as a cosmic hadron matter 
laboratory in Tables 1 and 2 . In each table the first column refers to the results of 
microscopic calculations which have the potential observational consequences spelled 
out in the second column, while the third column describes the current status of relevant 
observations. There is not time to discuss all aspects of these tables. In what follows I 
will focus primarily on the behavior of the crustal superfluid (for more details, and 
references, see Alpar et al., 1984a, and Pines and Alpar, 1985), which, as noted in 
Table 1, is believed to be responsible for 

• glitches, observed sudden changes in pulsar rotation rates, Qc and 6 C . 

• the long time scale (days to years) relaxation observed in pulsars following a 
glitch 

• the heating of old or cold neutron stars 

I shall then discuss briefly the physics issues raised in Table 2. 
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Table 1. Neutron Stars as Cosmic Hadron Matter Laboratories: 
The inner crust (p o /500 S pS po/2) 

Theoretical Consequences Possible observational Observations 

of hadron interaction consequences of hadron 

interaction 


• ‘'So crustal neutron 
superfluid (10 6 K5 T c 
10i OK), with 5* Rim, the 
size of a nucleus. 


• Pulsar spindown 
produces catastrophic 

vortex unpinning, 
leading to glitches 
(sudden pulsar spinups) 


• 16 giant glitches (10* 9 
5 AP/P 5 4x1 O' 6 ) seen 
in 7 pulsars. 


• Pinning of vortices in 
the crustal superfluid 
to the crust leads to lag 
of crustal superfluid in 
spinning down pulsars. 


• Following a glitch, 
vortex creep leads to a 
long time scale (d5t5y) 
recovery of Qc.“c. with 
(Ip/I) 5 10-2. 


• Long time scale 
recovery with (2d5t52y) 
seen for 12 glitches in 4 
pulsars. Excellent fit to 
vortex creep theory, with 
4 x 10-3 < Ip/I < 10-2. 


• Vortex creep also heats • Ti 929 + io rules out large 

old or cold neutron energy gaps and 

stars; observation of appreciable strong pinning 

pulsar surface tern- layers, 

peratures provides a 

direct measure of • Tcrab consistent with 

pinning strengths. standard cooling. 

• Direct determination of • Tyeia consistent with 

the internal pulsar unconventional cooling 

temperature is made plus vortex creep 

possi bl e by fits of post- h eati n g . 

glitch behavior using 

vortex creep theory. • T 0525-21 can be 

determined self- 
consistently from 
internal heating by 
vortex creep and fit 
to post-glitch behavior. 


• Whether the pinning is 
strong or weak depends 
sensitively on the mag- 
nitude of the pairing 
gaps. 
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Table 2. Neutron Stars as Cosmic Hadron Matter Laboratories: 
The quantum liquid core (p 0 ~ p ~ 3p 0 ) 

Theoretical consequences Possible observational Observations 

of hadron interaction consequences of hadron 

interaction 


• E.O.S. of neutron 
matter determines 
M-R relation, crustal 
extent, Mmax. Pmin (M)> 

Pmax(M) 


• 3 P 2 neutron superfluid 

• 1 S 0 proton superfluid 


• Pion condensates? 

• Quark liquids? 


Stellar wobble 

• 35 d cycle of Her X-1 

Limit on initial period 
of neutron stars, and on 
minimum period of 
millisec pulsars 

• Pms ~ 1 -5 ms? 

* P-I987A = 0.5ms? 

Starquakes? 

• Crab pulsar glitches? 

Prompt supernovae? 

• SN 1987A? 

Strong coupling 
between two super- 
fluids implies rapid 
crust-core coupling 

• No evidence for superfluid 
core lag in post-glitch 
behavior of pulsars 

Rapid cooling of stars 
with gapless pion 
condensates, quark 
cores. 

• Is the Vela pulsar too cold to be 
consistent with conventional 
theory? 


• Strange matter? 


Strange stars? 


Ruled out as "universal" by glitch 
observations. 
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The Crustal Neutron Superfluid 

The physical origin of the superfluidity of neutron matter is a BCS pairing instability 
brought about by the attractive hadron interaction at distances £ 1 fm. The superfluid 
transition temperature is given by 

T c ~ E f exp - 1/N(o)V ef f (1) 

where Ef is the neutron Fermi energy, N(o) the density of states per unit energy, and V e ff the 
strength of the most attractive pairing interaction. For densities p 5 p 0 , one does not sample 
appreciably the repulsive core (r c 5 1 fm) of the hadron interaction, and the pairs which make 
up the condensate will be in a state. At higher densities the tensor force provides the 
major source of attraction, and the pairs condense in a 3 P 2 state. Quite generally T c may be 
expected to be maximum at densities p ~ po/10. At lower densities, the density of states is too 
low to enable neutrons to take maximum advantage of the attractive hadron interaction, while 
at higher densities the neutrons begin to sample the repulsive core, so that V e ff falls off more 
rapidly with increasing density than N(o) increases. In Figure 2 I compare two calculations of 
the transition temperature; the upper curve represents early results obtained using separable 
potential or variational calculations; the two lower curves bracket some very recent results 
obtained in Urbana with a many-body theory which takes into account both short-range 
correlations produced by the repulsive core and long-range correlations produced by the 
induced particle-hole interactions (Ainsworth, Pines, and Wambach, 1989, and references 
therein). One test of our ability to use neutron stars as a hadron physics laboratory is the 
extent to which observations enable one to decide between these different theoretical results. 
This test is met in part: it is now possible to use observational results on the upper limit of the 
surface temperature of a relatively old pulsar to rule out the high transition temperatures 
characteristic of variational calculations. 

For pulsars rotating at an angular velocity Cl, it is energetically favorable for the 
neutron superfluid condensate to develop a vortex structure which enables it to mimic the 
effects of solid body rotation. One thus has lines of vortices parallel to the axis of rotation, 
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around which circulation is quantized in units of h/2m, of areal density such that the 
superfluid velocity, v s = £2r. The vortex density n is therefore n(h/2m) = 2Q. The size of the 
core of a vortex line is the superfluid coherence length 

p = 2E F , 13k F (fm-i fm 

7 tAk F A(Mev) (2) 

where kf is the Fermi wave vector and A the pairing energy gap; ^ is small compared to the 
average distance between vortices, 

d~3x10- 2 Q 1/2 cm. (3) 

In the stellar crust, vortices in the neutron superfluid coexist with a lattice of neutron- 
rich nuclei, of size Rn (~7 fm), spaced a distance b z apart. Since the size of the vortex core, 
is comparable to the size of a nuclear cluster (2Rn ~ 14 fm), it is energetically favorable for a 
vortex line to pass through a nuclear cluster if the energy gap for superfluid neutrons inside 
the cluster is smaller than outside, or to avoid the cluster if the gap for the cluster neutron 
superfluid is greater than the gap for the exterior crustal neutron superfluid. In the former 
case, it costs less energy to create a vortex line if it passes through the nuclear cluster; in the 
latter case, it would cost more. We can speak of the pinning (or threading) of vortex lines in 
the stellar crust; since the density of the neutron superfluid inside clusters is greater than that 
outside, we expect to encounter threading (avoidance of clusters) for crustal superfluid 
densities less than the density at which A(p) is maximum (~po/10) and pinning at higher 
densities. The pinning is intrinsic , in that the vortex core does not need a crystalline 
imperfection to pin, as is the case for a flux line in a typical type II superconductor. 

A rotating superfluid can change its angular velocity, Q, only through the motion of 
vortex lines; pinning therefore acts to prevent the superfluid from following the changes in the 
crustal angular velocity, Q c . produced by pulsar spin-down. As a result, a lag, co s £2 - £2 C , 
develops between the neutron superlufid and the crust as a pulsar spins down. This lag 
produces in turn a force, the Magnus force, on the vortices at density p and position r, 


fM = (Q s - Qc) P Kr 5 ^picr, ( 4 ) 

where k = h/2rriN is the quantum of vorticity. As the pulsar spins down, the Magnus force 
increases until eventually, at some critical angular velocity lag, o) cr , the vortices unpin and 
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move outwards. Because the distribution of pinning sites is not uniform, after traversing 
distance 8r, the vortices will repin, and the entire process begins anew, either in that, or 
another, portion of the stellar crust. This, then, is the basic glitch scenario: the storage of 
rotational energy in the superfluid as a result of pinning, followed by its sudden release: with 
the process repeating itself as the pulsar spins down, and vortices move in to replace those 
which have unpinned. 

The critical angular velocity lag is determined by balancing the Magnus force against 
the force per unit length of vortex line which must be overcome for vortices to move, 


f maxi.t Ep(A,p) 

' P (P)= WA,p)' 


(5) 


where E p is the pinning energy per nuclear cluster and b is the average spacing between 
pinning centers along the vortex line. Thus one finds 


cocr = 


picrb£ 


where the pinning energy per cluster, E p , is given by 

3 A 2 (p) 

E > (A - p) -tw^r 


( 6 ) 

(7) 


and y >s a factor of order unity which takes into account the overlap of the vortex core and the 
cluster, and the differences in condensation energies within and without the nuclear cluster. 

The pinning parameters are quite sensitive to the superfluid energy gap, A(p), for the 
density region of interest, 1/10 5 (p/p 0 ) 5 1/2. The physical quantities of interest in this region 
are the following: 


Nuclear spacing: 


30 fm 5 b z 5 50 fm 
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Neutron Fermi wavevector: 
Nuclear cluster size: 

Superfluid coherence length: 

Pinning energy: 

Energy to displace a nucleus a 
distance £ from equilibrium site 


0.82 fm-i 5 k f ~ 1 .2 frrv 1 
Rn ~ 7 fm 

% (fm) ~ 13 k f (fm-i)/A(MeV) 

E p (MeV) ~ k f (fm- 1 ) A 2 (MeV) 

72 e 2F2 0 

E L (MeV) ~ bk,(fm- 1 )/A 2 (MeV)(b z / 50 ) 3 - 

h ■ 


Depending then on A(r), pinning comes in one of three strengths. 


• Strong pinning 

It is energetically favorable to displace a nuclear cluster so that it pins more effectively 
to the vortex line. The condition for this is E p 2; El. In this regime, % s Rn, while the 
distance between pinning centers is b s ~ b z , from which it follows that the criterion for 
being in the strong pinning regime is 

A s (p) 2: 1 .5 kf 1 /4(fm)(b z /50)-3/4 MeV (8) 


and 


co| r (p) 2; 6.5 k f (fm- 1 )/pi 4 r6A(MeV)(b z /50)4 


( 9 ) 


where pi 4 is the density in units of 10 14 gm cm- 3 , and r6 the radial distance in units of 
10 6 cm. 




Weak pinning 

Vortices pin to whatever nuclear clusters are encountered. The spacing between 
pinning centers is then 


bw = bf/7t£ 2 = 230 


(b 2 /50) 3 A 2 (MeV) fm 

kf 


( 10 ) 


and 


C 0 cr s 0.4 rad s- 1 

r 6 kf(fm- 1 )(b 2 /50fm) 3 

is typically an order of magnitude smaller than co| r (p).. 


( 11 ) 
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Weak pinning exists for coherence lengths \ such that 

~ b z /2 ; (12) 

in other words, a vortex pins to only one nuclear cluster at a time. 


• Superweak pinning 

When the coherence length becomes comparable to, or larger than, the average 
spacing between nuclear clusters, 

£SW> bz /2 (13) 


vortices pin weakly to two or more adjacent nuclear clusters. Under these 


circumstances, it is obviously easier for a vortex line to move; there is no simple way to 
estimate either the effective pinning energy, (Ep)^, per cluster, or the 

Oil/ 

corresponding critical lag, co ^ ; we may, however, expect that 


(Ep) 


SW «E W 
eff p 

SW w 

“cr <<tD cr 


(14) 

(15) 


Vortex Creep 

The postglitch behavior and heating of old neutron stars is associated with the 
thermally activated motion of vortex lines. When Q s = &c, a vortex line sees a mixture of 


potential wells (the nuclear clusters) and barriers (produced by the Bernouilli force); as a 
result, its average drift velocity in the radial direction, <v r > = 0. When, however, there is a 

lag between the crust and the superfluid, the Bernoulli barriers are preferentially reduced in 
one direction, so that vortices will creep at an average velocity 


<v r > = v„exp-|f^ (16) 

where v 0 is a characteristic approach velocity (~ 10 7 cms' 1 ). This vortex creep is analogous 
to the thermally assisted creep of flux lines in a hard (Type-ll) superconductor studied by P.W. 
Anderson and his collaborators, and now of very great interest for the high temperature 
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superconductors. (Indeed, as Pethick and I have suggested elsewhere (Pethick and Pines, 
1987), the short coherence lengths found in the latter materials make possible intrinsic 
pinning of the kind discussed here for neutron stars.) As the lag, co, approaches co cr , such 
vortex creep enables the pinned superfluid to flow with greater ease, and so keep up with the 
spin-down induced decrease in £2 C . 

Quite generally, then, the crust is subject to both an external (magnetospheric) torque, 
Next, which produces spin-down, and an internal torque produced by the lagging crustal 
superfluid. The crustal angular velocity obeys the equation of motion, 

Ic^c = Next " (17) 

where l c is the total effective inertial moment of the crust plus the quantum liquid core (since 
the latter is found to couple to the crust on time scales (< 100s) which are short compared to 
those typically of interest for pulsar period irregularities), l p is the moment of inertia of the 
pinned crustal superfluid, and « is determined by the radial drift of the vortex lines, 


2Q 

q = -— <v r > 


(18) 


On combining the above equations, one obtains an equation for the lag, co = Q - Q c , 

&=|^|{1-[sinh^£f]/r|} , (19) 


where the vortex creep response parameter, given by 


ti = 


_ Iftblr 


4QqV 0 


exp 


kT 


( 20 ) 


determines the character of the superfluid response. 

Following a fllilctt, the resulting behavior of (or a>) is determined by the glitch- 
induced creep of the vortex lines produced by the sudden jump in Q c . The nature of that 
superfluid response is in turn determined by t \. Where rj « 1 , one has co « coc r , and the 
superfluid responds linearly to the perturbation, 5co(o), produced by a glitch; one finds 

co(t) = cobo - Scofoje^i 


( 21 ) 
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where the linear response relaxation time, 

t. - kT ^cr exD 

’‘4E pncVo 9XP kT (22) 

depends exponentially on the pinning energy and stellar temperature. On the other hand, 
where rj » 1 , one has co ~ co C r, and the crustal superfluid response depends non-linearly on 
the initial conditions. One finds 

= — I ) 

lc 1 + [exp(^_) - 1] exp- t/xni (23) 

where t 0 is an offset time which reflects the non-linear dependence on initial conditions, 


_ Mo) 

'"ISP 


and the non-linear response time, 


t - kT tQcr _ kT 1 

E p |0.| pkrb^ (25) 

depends linearly on the temperature and is independent of the pinning energy. In the weak 
pinning regime, one finds 

Tf), = 55 days kT(k ^- v) - 3 \ — 

A(Mev) I^J^Veial r 6 (b z /50fm) 3 [k,(fm-i)] 2 ( 26 ) 

on normalizing the pulsar spindown rate to that of the Vela pulsar, O/eia = 10' 10 rad s -2 . 

The transition from linear to non-linear response occurs at T n ® = t*, or equivalently, at a 

transition ratio of the pinning energy to the stellar temperature, 


= = 35 + in (ts/1 OM (27) 

where ts j (W/2|£2|<x> is the pulsar spindown time. 

To the extent that one knows the pulsar internal temperature and the superfluid energy gap 
and pinning energy as a function of density, one can then determine regions of linear and 
non-linear response in a given pulsar. Conversely, if the post-glitch behavior of a given 
pulsar can be shown to be linear or non-linear, one can use that information to deduce the 
character of A(p) (Alpar, Cheng, and Pines, 1989). 
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Energy dissipation by vortex creep is significant when conditions are such that 
appreciable portions of the neutron star are in the non-linear creep regime, for which o> ~ co cr . 
Under these circumstances, the energy dissipation rate in steady state is 

Ediss = Ip©crl^cl (28) 

This dissipation is large enough to be observable in the surface luminosity of older pulsars 

which have radiated away much of their initial heat content. 

Comparison between theory and observation 

In considering the observational evidence for superfluidity in neutron stars, one can 
ask questions at several different levels. First, can one uniquely deduce the presence of 
pinned crustal neutron superfluid from observations of pulsar glitches (eight in the Vela 
pulsar, three in the Crab pulsar, one each in pulsars 0525+21, 0355+54, 1641-45, 1325-43 , 
and 2224+65) and postglitch response. Here my answer would be an unambigious, "Yes." 
The evidence that glitches originate in the interior of neutron stars is overwhelming, and 
pinning of the crustal neutron superfluid is the only viable candidate. The magnitude of the 
jumps in both £2 C and Q; is easily understood as resulting from the catastrophic unpinning of 
vorticity in the crustal neutron superfluid involving pinning regions with inertial moments, (Ip/I) 
5 10' 2 , in accord with our theoretical expectations. Although we lack as yet a detailed theory 
of how such catastrophic events are triggered, the presence of inhomogeneous pinning 
regions makes glitches inevitable, while the time between the frequent glitches observed for 
the Vela and Crab pulsars is in accord with that calculated for pulsar spin down to enable 
unpinned vortices in these pulsars to move into pinning regions in sufficient number to initiate 
a subsequent glitch. 

Vortex creep theory likewise provides a natural explanation for the long time scale 
recovery of Qc and £2 C following a glitch; the observed time scales for that recovery, which 
range from a few days to a few months, are those expected for linear or non-linear response 
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in weak or superweak pinning regions. The four most recent Vela glitches and the latest 
Crab pulsar glitch have all been caught within a day of the glitch occurrence; fits to the 
postglitch data can be made with combinations of simple exponential decays, with two 
distinct time scales of a few days and a few months respectively, combined with the persistent 
shift in |Gt c | which is the hallmark of a substantial region of non-linear vortex creep. 

Is it possible to proceed to a second level of comparison of theory and observation, 
and use vortex creep theory to draw from observation some quantitative conclusions 
concerning the magnitude of critical angular velocities and pinning energies? Again the 
answer is "yes". For example, the failure to detect the thermal luminosity of the radio pulsar 
PSR 1929+10 enabled Alpar et al. (1987) to place an upper limit on vortex creep heating of 
this pulsar, and to conclude that <co C r> ~ 0.7 rad s’ 1 ; hence there can be no appreciable 
regions of strong pinning in this neutron star (and by inference, in any neutron star.) Alpar, 
Cheng and Pines (1989) have recently reviewed the application of vortex creep theory to 
postglitch behavior and arrive at the following conclusions: 

• There must be present, in the Crab pulsar, substantial regions (Ip/I £ 10' 3 ) in which the 

w 

pinning energy, Ep £ 1 MeV; these are the weak pinning regions expected for energy 
gaps, A, such that 1 MeV 5 A 5 1.5 MeV. 

• There must be present in PSR 0355-44, for which a linear postglitch response with a 

relaxation time of ~ 44d has been observed, appreciable regions of superweak 

sw 

pinning, with pinning energies, E p £ 50 kev. 

• The observed relaxation time 5 1 00d for the Vela pulsar could reflect either the non- 
linear response of regions of superweak pinning or the linear response of regions of 
weak pinning. 

• For PSR 0525-21 , the long relaxation time, 1 50d, extracted from a fit to the postglitch 
data, reflects the non-linear response of a region of weak pinning and is consistent 
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with an internal temperature of ~ 1.3x10‘ 2 kev calculated assuming vortex creep 
heating of this old (£ 1.5x10 6 y) pulsar. 

Given the difficulty in carrying out accurate calculations of energy gaps in neutron 

matter, and ambiguities in fitting postglitch response with linear or non-linear regimes, the 

agreement between vortex creep theory based on the APW gap calculations and observation 

is surprisingly good. As of this writing, a reasonable working hypothesis is that one has, for 

w 

all neutron stars, a region of weak pinning, Ep ~ IMeV, extending from densities of order 

3x10 13 gm cm’ 3 to 9x10 13 gm cm -3 , followed (in density) by regions of superweak pinning, 
svv 

E p 5 50kev, extending from densities of order 9x10 13 gm cm' 3 to ~ 2x10 14 gm cnr 3 . 

The Quantum Liquid Core 

Among the principal physics issues connected to the quantum liquid core of neutron 
stars are the following: 

• The character of the equation of state of neutron matter at densities greater than p 0 , 
the density of nuclear matter. 

• The existence of neutron and proton superfluid states. 

• The possible existence of a pion condensate (for p £ 2p 0 ). 

• A possible transition from hadron matter to a quark liquid at p £ 3p 0 . 

• The possible existence of strange stars, made up of strange quark matter. 

The first issue, although it is not likely to have appealed to Landau (who would likely have 
regarded it as undecidable and hence uninteresting), is of considerable importance for 
understanding neutron stars, since the neutron matter equation of state at densities in 
excess of p 0 determines the mass-radius relation, the crustal extent, the maximum mass, the 
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minimum rotational period, and the maximum crustal density for a given stellar mass. As 
Pandharipande et al. (1976) showed some time ago, one can construct a variety of models 
for the neutron interaction which are consistent with terrestrial constraints, and which lead to 
equations of state (EOS) which span a range from "soft" EOS (derived from models for which 
the average neutron interaction energy is attractive at nuclear densities) to "stiff" EOS 
(derived from models for which the average neutron energy is repulsive at sub-nuclear 
densities). In contrast to the 1 .33Me star depicted in Fig. 1 , which has been calculated using 
a moderately stiff equation of state, a star of comparable mass calculated with a soft 
equation of state based on the phenomemological Reid model for neutron-neutron 
interaction will have a much smaller radius (~8km), a significantly smaller crustal extent (5 
800 meters), and higher central density (~ 2 x 10 15 g crrr 3 ). The minimum rotational periods 
for stars of ~ 1 .4 Mq also differ substantially, ranging from ~ 1 .3ms for a stiff EOS, to ~ 0.7ms 
for a Reid EOS. 

Observational evidence in favor of a stiff EOS comes from the identification by Trumper 
et al. (1 986) of the 35d cycle of the pulsating x-ray source. Her X-1 as originating in free 
precession of the rotating neutron star, analogous to the Chandler wobble of the earth, and 
from the minimum magnitude of the inertial moment of the pinned crustal superfluid required 
to explain postglitch behavior (see Pines, 1987, for a discussion.) Should the minimum 
rotational period of a millisec pulsar be ~ 1.5ms, as Friedman et al. (1988) and Lipunov and 
Postnov (1 988) have suggested, again a quite stiff EOS is required. Under these 
circumstances one would be unlikely to find neutron stars with quark cores, since the 
calculated central density of a maximum mass neutron star is still rather less than 3p 0 , the 
most plausible density at which the hadron matter-quark liquid transition is likely to occur. It 
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is also worth noting that such stars would possess a sufficiently substantial crust that, when 
young, they would be subject to comparatively frequent starquakes induced by pulsar 
spindown; it is possible that such quakes could be responsible for the Crab pulsar glitches 
(Pines, 1987). 

On the other hand, as discussed recently by Brown (1988), current calculations show 
that "prompt" supernovae in which the stellar envelope is blown off immediately after 
implosion of the core, are only possible if the EOS remains soft up to densities around 4p 0 ; 
thus a stiff EOS favors "delayed" supernovae, in which neutrons act for some time following 
collapse to drive off the stellar envelope. 

Should the 0.5ms periodicity reported for the expected neutron star remnant of 
SN1987a (Kristian et al. 1989) be confirmed and identified with the rotational period of the 
star, "almost all the main features of current theories of neutron star formation, structure, 
evolution, and, most fundamentally, current models of nuclear interactions and the equation 
of state of matter at the highest densities" will require drastic revision (Alpar et al., 1 989). As 
many authors (Haensel and Zdunik, 1989; Friedman et al., 1989; Ipser and Lindblom, 1989; 
Sata and Suzaki, 1989; Shapiro et al., 1989) have now demonstrated, the existence of a 
uniform rotating star with a rotational period less than 0.5ms places unusually severe 
constraints on the equation of state, requiring that it be remarkably soft and lie within a 
pathologically narrow range of neutron star models. Thus the question considered above, of 
"soft" vs "stiff" equations of state becomes moot; all of the standard models presently in use 
are ruled out, and in all likelihood the pulsar would collapse into a black hole upon 
spindown. 
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In conclusion, (and assuming that 0.5msec pulsars do not exist!) I should like to 
address briefly the other three issues mentioned above. Pion condensates in neutron stars 
at densities greater than ~2p 0 continue to be a viable theoretical possibility; the best 
observational handle we have on their existence is the anomalously rapid (compared to the 
presence of, say, superfluid neutron matter) cooling of a neutron star such condensates 
make possible if they are in part gapless. Do there exist pulsars which appear to be far 
colder than conventional cooling theory would predict? Thus far the most promising 
candidate for such a star is the Vela pulsar. Ogelman and Zimmerman (1988) have 
deduced an interior temperature of ~ 1 1keV from their recent EXOSAT observations of its 
surface temperature; such a temperature is somewhat lower than the predictions (~16keV) of 
most standard cooling theories. A still lower temperature (~ 1.4keV) was obtained by Alpar 
et al. (1984) who interpreted the observed 60d postglitch relaxation time as reflecting weak 
pinning and non-linear response. As noted earlier, should the postglitch response be linear, 
or reflect a superweak pinning regimes, there is no need to invoke such a low temperature. 

Witten's proposal that neutron stars might in fact be composed of strange quark matter 
was an intriguing one; however Alpar (1987) pointed out that since such stars would almost 
certainly not possess a large enough solid outer crust to produce the observed amplitudes of 
glitches and post-glitch relaxation, it can at least be concluded that not all pulsars are 
strange. Indeed, to the extent that glitch statistics are consistent with the hypothesis that all 
pulsars glitch occasionally, perhaps no pulsars are strange. 

Finally, what evidence do we have for the superfluidity of the neutron and proton 
liquids in the quantum liquid core? Because the core superfluid is now believed to couple on 
a very short time scale (~ min.) to the stellar crust (Alpar and Sauls, 1989), one has no direct 



20 


evidence of superfluidity from glitch observation. Since, however, Hamilton et al. (1989) 
were observing the Vela pulsar at the time of its eighth glitch, it is possible that the shortest 
time scales found in these observations will provide direct evidence concerning the existence 
of a superfluid neutron core, as well as providing invaluable observational limits on the 
duration of the glitch itself. Indeed, for "glitch theorists" this observation may well play a role 
comparable to that of SN1987A for "supernova theorists." 

It gives me great pleasure to thank Ali Alpar for stimulating discussions on these and 
related topics, and for a critical reading of this manuscript. 
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Figure 2. Calculations of the superfluid neutron transition temperature. The upper curve 
represents results obtained using a separable potential or variational approach; 
the two lower curves, APW-1 and APW-2, represent a bound on recent many-body 
calculations (Ainsworth, Pines, and Wambach, 1989) which take into account both 
short range correlations produced by the repulsive core and long-range 
correlations produced by induced particle hole interactions. 



